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cross  sectional  area,  m 
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1.0  EXECUTIVE  SUMMARY 


Often  it  is  desirable  to  cool  electronic  components  to  minimize  the  temperature  gradient 
across  the  heat  generating  surface.  This  isothennal  requirement  is  often  required  to  in  order  to 
minimize  thennal  stresses,  in  the  case  of  power  electronics,  or  in  the  case  of  laser  components, 
such  as  laser  diodes,  to  minimize  optical  distortion. 

The  objective  of  this  analysis  is  to  detennine  the  effects  of  boundary  condition(s)  that  will 
influence  the  steady-state  temperature  gradient  along  the  top  heated  surface  of  cylindrical  and 
rectangular  posts.  This  is  accomplished  by  casting  the  dimensionless  analytical  solutions  in 
graphical  form  similar  to  Heisler-Grober  Charts.  Both  dimensional  and  dimensionless 
analytical  solutions  are  presented  followed  by  graphical  presentations  in  dimensionless  form. 

In  this  manner,  the  “best”  boundary  conditions  can  be  used  to  specify  the  cold  plate 
requirements  to  optimally  cool  the  electronics  package.  The  generated  dimensionless  plots  also 
allow  for  easy  graphical  detennination  of  these  boundary  conditions  as  defined  in  a 
dimensionless  form  which  can  be  readily  interpreted  and  used  to  determine  optimal 
dimensional  geometries  and/or  thermo  physical  material  properties  to  minimize  temperature 
gradients  across  a  heated  electronic  surface. 
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2.0  INTRODUCTION 


Often  it  is  desirable  to  cool  electronic  components  in  such  a  manner  as  to  minimize  the 
temperature  gradient  across  the  heat  generating  surface.  This  isothermal  requirement  is  often 
required  to  in  order  to  minimize  thermal  stresses,  in  the  case  of  power  electronics,  or  in  the 
case  of  laser  components,  such  as  laser  diodes,  to  minimize  optical  distortion.  Traditionally, 
electronics  cooling  has  been  accomplished  by  mounting  cold  plates  or  cooling  fins  to  the 
electronics  package.  In  doing  so,  the  “best”  cooling  approach  was  typically  taken  from  a 
thermal  management  “tool  box”  consisting  of  proven  concepts.  This  bolt-on  approach  to 
electronics  thermal  management  is  only  partially  effective  and  does  little  to  optimize  the 
thermal  management  of  electronics  with  regard  to  minimizing  thennal  stresses  within  the 
electronics  package.  Another  option  to  electronics  thermal  management  is  to  define  the 
boundary  conditions  in  order  to  minimize  thennal  stresses  within  the  electronics  package 
generating  the  heat.  The  specific  cooling  approach  is  then  designed  to  optimally  match  the 
boundary  conditions  that  will  minimize  thennal  stresses  within  the  electronics  package.  Once 
this  is  done  one  can  account  for  the  necessary  thermal  management  components  needed  to 
transport  the  heat  to  the  final  heat  sink.  Yovanovich1  investigated  thennal  spreading  resistances 
in  a  finite  circular  cylinder  with  side  and  end  cooling.  He  developed  a  general  dimensionless 
analytical  solution,  for  a  circular  heat  source,  showing  that  the  dimensionless  thermal 
resistance  depended  upon  four  dimensionless  system  parameters;  relative  heat  source  size, 
relative  cylinder  thickness,  and  side  and  end  surface  Biot  numbers.  Muzychka  followed  up  on 
the  research  of  Yovanovich  and  examined  the  dimensionless  analytical  solutions  for  the 
thermal  spreading  resistance  of  circular  cylinder  and  rectangle  geometries.  These  solutions 
addressed  geometries  consisting  of  both  isotropic  and  compound  systems  and  resulted  in 
simple  expressions  to  detennine  the  impact  of  edge  cooling  on  the  thermal  spreading  resistance. 
However,  these  studies  do  not  take  into  account  temperature  variation  at  the  heated  surface. 

The  objective  of  this  analysis  is  to  detennine  the  effects  of  boundary  condition(s)  that  will 
influence  the  steady-state  temperature  gradient  along  the  top  heated  surface  of  cylindrical  and 
rectangular  posts.  This  is  accomplished  by  casting  the  dimensionless  analytical  solutions  in 
graphical  form  similar  to  Heisler-Grober  Charts.  '  Both  dimensional  and  dimensionless 
analytical  solutions  are  presented  followed  by  graphical  presentations  in  dimensionless  form. 


3.0  ANALYTICAL  FORMULATION 

A  simplified  example  of  this  analytic  approach,  in  rectangular  and  cylindrical,  dimensional 
and  dimensionless  forms,  is  shown  in  Figure  1 .  In  this  example,  a  uniform  and  constant  heat 
flux  is  applied  to  the  top  surface  of  a  solid  material.  The  centerline  of  the  material  lies  along  y 
=  7=  0orr  =  f?  =  0  and  is  defined  as  an  adiabatic  surface.  A  convective  boundary  of  the  third 
kind  lies  along  y,r  =  b  or  Y,R  =  1  and  is  bounded  by  temperature,  Tm.  The  base  of  the  material 
represents  a  cold  plate,  that  is  a  boundary  condition  of  the  first  kind,  at  a  fixed  constant 
temperature  T  =  Tc  or  6  =  1 . 
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dT  Q!  A 
dz  k 


-  Ill  HI 

dy 


(a) 


T  =  TC 


dO 

dZ 


=  -ri 


(b) 


0  =  1 


dO 

dY 


+  Bi0  =  0 


dT  Q/A 
dz  k 


d T 
dr 


dT 

dr 


+  HT  =  HT 


(C) 


T  =  T 


d0 

dZ 


=  ~0 


—  +  Bi0  =  O 
dR 


Figure  1.  Example  for  (a)  rectangular,  dimensional  case;  (b)  rectangular,  dimensionless;  (c)  cylindrical, 
dimensional  case;  and  (d)  cylindrical,  dimensionless  case. 
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3.1  Rectangular 


The  dimensionless  fonn  of  the  example  shown  in  Figure  lb  is  formulated  using  the 
dimensionless  parameters 


(r-rj 

(L-tJ 


y 

and  Y  =  —. 
b 


The  dimensionless,  steady-state,  fonn  of  the  conservation  of  energy  now  becomes 

^4  +  ^4  =  0,  0  <Z<£,  0  <  7  <  1,  (la) 

dr2  dz2  * 

with  the  boundary  conditions 


—  =  0,  7  =  0, 
dY 

(lb) 

—  +  BiO  =  0,  7  =  1, 
dY 

(lc) 

dO  7  n 

—  =  -ri,  z  =  o, 

dZ 

(Id) 

and  0  =  1,  Z  = 

(le) 

The  dimensionless  solution  of  Eqs.  (la-le)  is  given  by 

-  (  y 

cosh(/?mZ)  +  ^sinh(/?m(£-Z))  - 

m=i  V  Pm 

where  pms  are  the  positive  roots  of 


0(Z,Y)  =  £ 


2cos(j3nY)(sni/Jm)(]3m2+Bi2)' 
j3m[(j3m2+Bi2)  +  Bi]Cosh(/3J)y 


,  (2) 


Pm  tan(AJ  =  BL 


3.2  Cylindrical 


The  dimensionless,  steady-state,  form  of  the  example  shown  in  Figure  Id  is  formulated 
using  the  dimensionless  parameters 


{t-tJ 

(K-tJ 
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r 

and  R  =  — . 
b 


(3a) 

(3b) 

(3c) 

(3d) 

(3e) 


Equations  2  and  4  were  solved  to  detennine  the  effect  on  the  surface  temperature  across  the 
top  surface  at  Z  =  0  for  both  the  rectangular  and  cylindrical  cases.  The  top  surface  temperatures 
at  the  centerline,  0(0, 0),  and  temperature  variation,  0(0,0)  -  0(0,1),  were  calculated  for  varying 
dimensionless  heat  input, ;/,  aspect  ratio,  <f,  and  Biot  number  of  the  fonn  1/0/.  The  eigenvalues 
were  solved  using  a  bisection  method.  The  number  of  eigenvalues,  used  in  the  series  expansion, 
was  chosen  in  such  a  manner  that  when  doubled,  the  resulting  solution  varied  less  than  1%. 

The  results  of  these  calculations  are  shown,  as  plots,  in  Figs.  2-15  and  Figs.  16-29  for  the 
rectangular  and  cylindrical  cases  respectively. 

There  were  small  variations  when  the  results  for  the  rectangular  and  cylindrical  cases  were 
compared;  but  overall,  both  show  the  same  trend  as  dimensionless  heat  input,  rj,  aspect  ratio,  <f, 
and  Biot  number  of  the  fonn,  1/Bi,  were  varied. 
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4.1  Rectangular 

Figures  2-15  are  grouped  into  three  sets  of  plots  each  highlighting  the  effects  of  the  three 
varying  parameters;  dimensionless  heat  input,  rj,  aspect  ratio,  2,  and  Biot  number  of  the  form 
1/Bi.  Figures  2-5  show  the  variation  in  the  top  surface  dimensionless  centerline  temperature, 
0(0,0),  and  temperature  variation,  0(0,0)  -  0(0,1),  plotted  against  1/Bi  with  aspect  ratio,  <f, 
varying  from  0.01  to  5  and  dimensionless  heat  input,  //,  varying  from  -0.5  to  0.5.  Figures  6-10 
show  the  variation  in  the  top  surface  dimensionless  centerline  temperature,  0(0,0),  and 
temperature  variation,  0(0,0)  -  0(0,1),  plotted  against  dimensionless  heat  input, ;/,  with  1/2?/ 
varying  from  0.01  to  1000  and  aspect  ratio,  varying  from  0.25  to  5.  Figures  11-15  show  the 
variation  in  the  top  surface  dimensionless  centerline  temperature,  0(0,0),  and  temperature 
variation,  0(0,0)  -  0(0,1),  plotted  against  aspect  ratio,  with  1/Bi  varying  from  0.01  to  1000 
and  dimensionless  heat  input,  //,  varying  from  -0.5  to  0.5. 


0.01  0.1  1  10  100  100(  0.01  0.1  1  10  100  1000 
1/Bi  1/Bi 

Figure  2.  Top  surface  dimensionless  centerline  temperature,  6>(0,0),  and  temperature  variation,  0(0,0)  -  0(0,1), 
plotted  against  1/0/  for  varying  aspect  ratio,  <f;  (a),  (b)  t]  =  -0.5,  (c),  (d)  //  =  -0.2,  and  (e),  (f)  //  =  -0.1. 
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(9(0,0)  (9(0,0) 


Figure  3.  Top  surface  dimensionless  centerline  temperature,  (9(0.0).  and  temperature  variation,  0(0,0)  -  0(0,1), 
plotted  against  1  /Bi  for  varying  aspect  ratio,  <f;  (a),  (b)  rj  =  0.1,  (c),  (d)  tj  =  0.2,  and  (e),  (f)  ij  =  0.5. 
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Figure  4.  Top  surface  dimensionless  centerline  temperature,  0(0,0),  and  temperature  variation,  0(0,0) 
-  0(0,1),  plotted  against  1/5/  for  varying  dimensionless  heat  input,  //;  (a),  (b)  =  0.25,  (c),  (d)  =  0.5, 
and  (e),  (f)  C=  1. 
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0(0,0) 


0.01  0.1  1  10  100  1000  0.01  0.1  1  10  100  1' 
1/Bi  J/Bi 


Figure  5.  Top  surface  dimensionless  centerline  temperature,  0(0,0).  and  temperature  variation,  0(0,0)  - 
0(0,1),  plotted  against  1  /Bi  for  varying  dimensionless  heat  input,  //;  (a),  (b)  C  =  2,  (c),  (d)  C  =  3,  (e),  (f)  C 
=  4,  and  (g),  (h)£  =  5. 
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Figure  6.  Top  surface  dimensionless  centerline  temperature,  0(0,0),  and  temperature  variation,  0(0,0)  - 
0(0,1),  plotted  against  dimensionless  heat  input,  //,  for  varying  1/0/;  (a),  (b),  (c),  (d)  £  =  0.25  and  (e),  (f), 
(g)>  (h)  C  =  0.5. 
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0(0,0)  0(0,0) 


Figure  8.  Top  surface  dimensionless  centerline  temperature,  0(0,0),  and  temperature  variation,  0(0,0)  - 
0(0,1),  plotted  against  dimensionless  heat  input,  //,  for  varying  1  tBi\  (a),  (b),  (c),  (d)  C  =  3  and  (e),  (f),  (g), 
(h)  c-  5. 
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Figure  9.  Top  surface  dimensionless  centerline  temperature,  0(0,0),  and  temperature  variation,  0(0,0)  - 
0(0,1),  plotted  against  dimensionless  heat  input,  //,  for  varying  aspect  ratio,  <f;  (a),  (b)  1/Bi  =  1000,  (c), 
(d)  1/Bi  =  100,  and  (e),  (f)  1/Bi  =  10. 
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Figure  10.  Top  surface  dimensionless  centerline  temperature,  (9(0,0),  and  temperature  variation,  (9(0,0) 
0(0,1),  plotted  against  dimensionless  heat  input,  //,  for  varying  aspect  ratio,  Q  (a),  (b)  1/Bi  =  5,  (c),  (d) 
1/Bi  =  1,  and  (e),  (f)  1/Bi  =  0.01. 
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61(0,0) 


f 


Figure  11.  Top  surface  dimensionless  centerline  temperature,  (9(0,0),  and  temperature  variation,  (9(0,0) 
6>(0,1),  plotted  against  aspect  ratio,  <f,  for  varying  1  /Bi;  (a),  (b)  rj  =  -0.5,  (c),  (d)  tj  =  -0.2,  and  (e),  (f)  tj  = 
0.1. 
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Figure  13.  Top  surface  dimensionless  centerline  temperature,  0(0, 0),  and  temperature  variation,  0(0,0) 
0(0,1),  plotted  against  aspect  ratio,  £  for  varying  dimensionless  heat  input,  rj;  (a),  (b),  (c),  (d)  1/Bi  = 
1000  and  (e),  (f),  (g),  (h)  1/Bi  =  100. 
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4.2  Cylindrical 


Similar  to  the  rectangular  coordinate  case,  Figs.  16-29  are  grouped  into  three  sets  of  plots 
each  highlighting  the  effects  of  the  three  varying  parameters;  dimensionless  heat  input, //, 
aspect  ratio,  and  Biot  number  of  the  form  1/2?/.  Figures  16-19  show  the  variation  in  the  top 
surface  dimensionless  centerline  temperature,  0(0,0),  and  temperature  variation,  0(0,0)  -  0(0,1), 
plotted  against  1/2?/  with  aspect  ratio,  <f,  varying  from  0.01  to  5  and  dimensionless  heat  input,  rj, 
varying  from  -0.5  to  0.5.  Figures  20-24  show  the  variation  in  the  top  surface  dimensionless 
centerline  temperature,  0(0,0),  and  temperature  variation,  0(0,0)  -  0(0,1),  plotted  against 
dimensionless  heat  input, ;/,  with  l/Bi  varying  from  0.01  to  1000  and  aspect  ratio,  2,  varying 
from  0.25  to  5.  Figures  25-29  show  the  variation  in  the  top  surface  dimensionless  centerline 
temperature,  0(0,0),  and  temperature  variation,  0(0,0)  -  0(0,1),  plotted  against  aspect  ratio,  2, 
with  1/2?/  varying  from  0.01  to  1000  and  dimensionless  heat  input,  rj,  varying  from  -0.5  to  0.5. 


3.5 


-1.5  3 - 1 - 1 - 1 - 1 - 1  -0.4 


Figure  16.  Top  surface  dimensionless  centerline  temperature,  0(0,0),  and  temperature  variation,  0(0,0)  - 
0(0,1),  plotted  against  l/Bi  for  varying  aspect  ratio,  <f;  (a),  (b)  t]  =  -0.5,  (c),  (d)  //  =  -0.2,  and  (e),  (f) ;/  =  - 
0.1. 
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Figure  17.  Top  surface  dimensionless  centerline  temperature,  0(0, 0),  and  temperature  variation,  0(0,0)  - 
0(0,1),  plotted  against  1  IBi  for  varying  aspect  ratio,  <f;  (a),  (b)  t)  =  0.1,  (c),  (d)  t]  =  0.2,  and  (e),  (f)  tj  =  0.5 
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Figure  18.  Top  surface  dimensionless  centerline  temperature,  0(0,0),  and  temperature  variation,  0(0,0)  - 
0(0,1),  plotted  against  1/Bi  for  varying  dimensionless  heat  input,  t]\  (a),  (b)  =  0.25,  (c),  (d)  =  0.5,  and 

(e),(f)C=L 
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Figure  20.  Top  surface  dimensionless  centerline  temperature,  0(0,0),  and  temperature  variation,  0(0,0)  - 
0(0,1),  plotted  against  dimensionless  heat  input,  rj,  for  varying  1/Bi;  (a),  (b),  (c),  (d)  £=  0.25  and  (e),  (f), 
(g)>  (h)  <f=  0.5. 
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Figure  21.  Top  surface  dimensionless  centerline  temperature,  0(0, 0),  and  temperature  variation,  0(0,0)  - 
0(0,1),  plotted  against  dimensionless  heat  input, for  varying  1  /Bi;  (a),  (b),  (c),  (d)  f  =  1  and  (e),  (f),  (g), 
(h)  C  =  2. 
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Figure  22.  Top  surface  dimensionless  centerline  temperature,  0(0,0),  and  temperature  variation,  0(0,0)  - 
0(0,1),  plotted  against  dimensionless  heat  input,  rj,  for  varying  1/5/;  (a),  (b),  (c),  (d)  <f=  3  and  (e),  (f),  (g), 
(h)C  =  5. 
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Figure  23.  Top  surface  dimensionless  centerline  temperature,  0(0,0),  and  temperature  variation,  0(0,0)  - 
0(0,1),  plotted  against  dimensionless  heat  input,  tj,  for  varying  aspect  ratio,  <f;  (a),  (b)  1/0/  =  1000,  (c), 

(d)  1/Bi  =  100,  and  (e),  (f)  1/Bi  =10. 
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Figure  24.  Top  surface  dimensionless  centerline  temperature,  0(0,0),  and  temperature  variation,  0(0,0)  - 
0(0,1),  plotted  against  dimensionless  heat  input,  //,  for  varying  aspect  ratio,  Q  (a),  (b)  1/Bi  =  5,  (c),  (d) 
1/Bi  =  1,  and  (e),  (f)  1/Bi  =  0.01. 
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Figure  25.  Top  surface  dimensionless  centerline  temperature,  0(0,0),  and  temperature  variation,  0(0,0) 
0(0,1),  plotted  against  aspect  ratio,  <f,  for  varying  1/Bi;  (a),  (b)  rj  =  -0.5,  (c),  (d)  ij  =  -0.2,  and  (e),  (f)  tj  = 
0.1. 
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4.3  Comparison  of  Rectangular  and  Cylindrical  Results 

Both  rectangular  and  cylindrical  results  compare  closely  with  only  slight  differences.  One 
aspect  noted  is  the  fact  that  there  are  some  unique  combinations  of  dimensionless  heat  input 
and  aspect  ratio  combinations  that  are  independent  of  l/Bi.  Figures  30  and  31  show  these 
dimensionless  heat  input  and  aspect  ratio  combinations  that  are  independent  of  l/Bi  for  the  top 
surface  dimensionless  centerline  temperature,  0(0,0)  and  the  top  surface  dimensionless 
temperature  variation,  0(0,0)  -  0(0,1).  These  combinations  can  be  readily  determined  from  Figs. 
6-8  and  Figs.  20-22,  for  the  rectangular  and  cylindrical  geometries  respectively,  as  the  location 
of  the  intersection  point  for  the  constant  l/Bi  lines.  In  the  plots  shown  in  Figs.  30  and  3 1  were 
generated  with  constant  lines  l/Bi  =  0.01  and  1000. 


Figure  30.  Rectangular  dimensionless  heat  input  and  aspect  ratio  combinations  that  are  independent  of 
l/Bi  for  (a)  top  surface  dimensionless  centerline  temperature,  0(0,0)  and  (b)  top  surface  dimensionless 
temperature  variation,  0(0,0)  -  0(0,1). 


(a)  -  (b) 


Figure  31.  Cylindrical  dimensionless  heat  input  and  aspect  ratio  combinations  that  are  independent  of  l/Bi 
for  (a)  top  surface  dimensionless  centerline  temperature,  0(0,0)  and  (b)  top  surface  dimensionless 
temperature  variation,  0(0,0)  -  0(0,1). 
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5.0  CONCLUSION 


This  analytical  study  presents  results  that  can  be  used  to  determine  the  effects  of  boundary 
condition(s)  that  will  influence  the  steady-state  temperature  gradient  along  the  top  heated 
surface  of  cylindrical  and  rectangular  posts  with  sidewall  cooling.  Both  dimensional  and 
dimensionless  analytical  solutions  derived  and  the  results  are  presented  in  dimensionless 
graphical  form  similar  to  Heisler-Grober  Charts. 

This  specific  example  demonstrates  how  it  is  possible  to  detennine  the  cold  plate  surface 
and  environmental  boundary  conditions  that  are  required  to  minimize  either  the  centerline 
temperature  or  temperature  variation  across  the  heated  surface.  In  this  manner,  the  “best” 
boundary  conditions  can  be  used  to  specify  the  cold  plate  requirements  to  optimally  cool  the 
electronics  package.  The  generated  dimensionless  plots  also  allow  for  easy  graphical 
determination  of  these  boundary  conditions  as  defined  in  a  dimensionless  form  which  can  be 
readily  interpreted  and  used  to  detennine  optimal  dimensional  geometries  and/or  thenno 
physical  material  properties  to  minimize  temperature  gradients  across  a  heated  electronic 
surface. 
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